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Abstract 

The spin-fermion model has long been used to describe the quantum-critical behav- 
ior of 2d electron systems near an antiferromagnetic (AFM) instability. Recently, 
the standard procedure to integrate out the fermions to obtain an effective action 
for spin waves has been questioned in the clean case. We show that in the presence 
of disorder, the single fermion loops display two crossover scales: upon lowering the 
energy, the singularities of the clean fermionic loops are first cut off, but below a 
second scale new singularities arise that lead again to marginal scaling. In addi- 
tion, impurity lines between different fermion loops generate new relevant couplings 
which dominate at low energies. We outline a non- linear cr model formulation of the 
single-loop problem, which allows to control the higher singularities and provides 
an effective model in terms of low-energy diffusive as well as spin modes. 
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1 Introduction 



The spin-fermion model is a low-energy effective model describing the inter- 
action of conductance electrons (fermions) with spin waves (bosons). It has 
been used, e.g., to describe the quantum critical behavior of an electron system 
near an antiferromagnetic instability [TI2I3] . An important example where this 
might be realized experimentally is in itinerant heavy-fermion materials [1]. 
By integrating out the fermions completely, a purely bosonic effective action 
for spin waves is obtained. This action is written in terms of a bare spin propa- 
gator and bare bosonic vertices with any even number of spin lines. The value 
of each bosonic vertex is given by a fermionic loop with spin-vertex insertions: 
in general, these are complicated functions of all external bosonic frequencies 
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and momenta. Hertz [T] and Millis [2] considered only the static limit of these 
vertices, i.e., setting all frequencies to zero at finite momenta. In this limit the 
4-point vertex and all higher vertices vanish for a linear dispersion relation, 
while they are constants proportional to a power of the inverse bandwidth if 
the band curvature is taken into account. For the AFM the dynamic critical 
exponent is z = 2 due to Landau damping of spin modes by particle-hole pairs. 
The scaling of these vertices in d = 2 under an RG flow toward low energy 
scales is marginal for the 4-point vertex while all higher vertices are irrelevant 
{d + z = 4 is the upper-critical dimension). Thus, a well-defined bosonic action 
with only quadratic and quartic parts in the spin field is obtained. 

Recently, Lercher and Wheatley [5] as well as Abanov et al. [6] considered 
for the 2d case not only the static limit of the 4-point vertex but the full 
functional dependence on frequencies and momenta. Surprisingly, they found 
that in the dynamic limit, setting the momenta to zero at finite frequency, the 
4-point vertex is strongly divergent as the external frequencies tend to zero, 
implying an effective spin interaction nonlocal in time. The higher bosonic 
vertices display an even stronger singularity [7]. 

To assess the relevance of the singular bosonic vertices, Ref. [7] considered 
the scaling limit uo with z = 2. In this limit, the bosonic vertices are 

less singular than in the dynamic limit but the related couplings are still 
marginal, i.e., they cannot be neglected in the effective bosonic action, in ap- 
parent contradiction to Hertz and Millis. Employing an expansion in a large 
number of hot spots or fermion flavors, Ref. [B] argues that vertex correc- 
tions are resummed to yield a spin propagator with an anomalous dimension 
rj = 2/N = 1/4 (for A^ = 8). At the same time, z = 2 remains unchanged up 
to two-loop order, i.e., the frequency dimension is given by = 2(1 — ri/2). 

The infinite number of marginal vertices renders the purely bosonic theory 
difficult to use for perturbative calculations. A relevant question is whether the 
above difficulty persists upon the inclusion of a weak static disorder potential 
present in real materials. To address this issue we insert disorder corrections 
into single fermionic loops and find two different crossover scales: at frequencies 
uj ^ l/r (i.e., much larger than the impurity scattering rate) and momenta 
q ^ 1/i (with mean free path i = vpT, where vp is the Fermi velocity), the 
fermionic loops resemble the clean case, while below this scale the singularity is 
cut off by self-energy corrections and the loops saturate. However, at yet lower 
frequencies, a second crossover scale u ~ l/^rkpi), q ~ l/{i^/kpi) appears 
where the loops acquire a diffusive form due to impurity ladder corrections and 
the related couplings again scale marginally, as in the clean case. Therefore, 
in an intermediate energy range disorder regularizes the singular vertices and 
appears to restore Hertz and Millis theory, while ultimately at the lowest scales 
the disordered loops are as singular as the clean ones, albeit with a different 
functional form: the linear dispersion of the electrons is replaced by a diffusive 
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form. We outline a non-linear a model formulation of the disordered single- 
loop problem which allows us to identify all disorder corrections which exhibit 
the maximum singularity, and provides an action for spin modes coupled to 
low-energy diffusive electronic modes, instead of the original electrons. 

Finally, while all disorder corrections to a single fermion loop lead to couplings 
which scale at most marginally, impurity lines connecting different fermion 
loops are a relevant perturbation in = 2. We find that these diagrams may 
dominate the single-loop contributions below a; ~ 1/r, depending on the typ- 
ical values of the bosonic momenta. 

We proceed as follows: in the remaining part of this section, we introduce the 
model and the scaling arguments for the clean case. We then insert disorder 
corrections into a single fermion loop and discuss a class of most singular dia- 
grams in section [21 Their scaling behavior and the emergence of two crossover 
scales is the subject of section El The multi-loop diagrams are discussed in 
section HI Appendix |X] contains the non-linear a model for the disordered 
single-loop case. 

1.1 The spin-fermion model in 2d 

The 2d spin-fermion model is defined by the action 

S[^, i;, 0] = {^P, Go V) + (0, Xo V) + 9<Pi^^ 

for a fermionic field ^p, ijj and a bosonic spin field (/.Q The inverse fermionic 
propagator is Gq ^(ie, p) = ie — ^p in terms of the Matsubara frequency ie and 
a dispersion relation with a roughly circular Fermi surface (FS), which we 
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approximate by a quadratic dispersion = ^ /i with electron mass rrie, 

chemical potential n, Fermi momentum kp = y/2mefi, and constant density 
of states 27rpo = ^p/vp. Xo{q) is the bare spin propagator. 

We assume that the above model describes an AFM quantum critical point 
at finite q = Qc- The Fermi surface has so-called hot regions connected by 
exchange of g^, and cold regions where scattering off spin waves is weak. Here 
we shall assume an underlying lattice and a commensurate = (tt, tt), which 
is equivalent (up to a reciprocal lattice vector) to — q^,- 

When computing fermionic loops with only spin-vertex insertions, the momen- 
tum integration can be reduced to the region around two hot spots separated 
by q^, which we shall denote by a and a (Fig. [T|). The fermionic dispersion 

^ The explicit spin structure of the spin-fermion vertex is not relevant for the fer- 
mionic loops and needs to be specified only when spin lines are contracted. 
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Fig. 1. The Fermi surface with the hot spots separated by the wave vector q^. 
relation is linearized around any hot spot a at momentum pg'g as [6] 

where p denotes the distance from the hot spot. The components {vy) of 
the Fermi velocity vp parallel (perpendicular) to at a given hot spot a are 
related by ff^ = v'^ + Vy and Vx/vy = tan(0o/2), with 0o the angle between hot 
spots a and a as seen from the center of the circular Fermi surface. The case 
(1)0 = n {vy = 0) corresponds to perfect nesting, but here we consider a generic 
00 without nesting. For a pair of hot spots, the momentum integration can be 
written as 

where and are two independent momentum directions at hot spot a 
coincides with the radial direction at hot spot a), = Air'^Vp sin (po is the 
corresponding Jacobian which depends on the shape of the Fermi surface and 
the filling, and one still has to perform the summation over all = 8 hot 
spots. 

1.2 Clean fermionic 2n-loops 

The fermionic loops with 2n spin- vertex insertions — i.e., the 2?T,-point functions — 
are in general complicated functions of the external frequencies and momenta. 
The loop with two spin insertions contributes to the self-energy for the spin 
propagator and has the well-known Landau damping form for small frequen- 
cies u and momenta near [8|, 

J:{iuj,q^ qj = -7 \uj\ , 
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with the dimensionless strength of the spin fluctuations |6] 



7 = 271 JNg^' 



27cvp sin I 



4:71 V^Vy 



The inverse spin propagator resummed in the random-phase approximation 
has dynamical exponent z = 2, 



X {iiu,q) = m + 'y\uj\ + iy\q 



|2 

Qc\ , 



where the mass term m measures the distance from the quantum critical 
point and u ~ g'^/tp- Near criticality m 0, the momentum exchanged by 
scattering off a spin wave is peaked near q^. From here on, we shall denote 
by q the deviation from q^. In the commensurate case there are logarithmic 
singularities in the clean bosonic self-energy (from contracting two spin lines 
diagonally in Fig. [2]) that may lead to an anomalous dimension of the spin 
propagator [B]. 




— luji, q^ 



iu + iui, qc + q + qi 



-lUJ — IU2^ 



Q2 



Fig. 2. The clean fermionic four-loop 64, indicating the notation for the external 
frequencies and momenta. 

The 4-point function is given by the fermion loop with four spin insertions, 
which provides the bare two-spin interaction: 



, \uJi + Uj\ + \uJi — u\ — \UJ2 + Uj\ — \UJ2 — Uj\ 

where we have labeled the three independent external frequencies and mo- 
menta as shown in Fig. [2j This is a nonanalytic function whose value for 
— >^ 0, qr — > depends on the order of the limits: it vanishes in the static 
limit {uj first) while it diverges as 1/u) in the dynamic limit (q — > 
first). There is an important difference from the forward-scattering loop with 
small external momenta [9]: here, symmetrization of the external lines does 
not lead to loop cancellation, i.e., a reduction of the leading singularity. In- 
stead, symmetrization only modifies the prefactors but does not change the 
scaling dimension. 
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1.3 Scaling of the clean fermionic loops 

We recall the scaling behavior of the 2n-point functions [7j. Since no cancel- 
lation of the leading singularity occurs, we only need to consider the power of 
external frequency and momentum and not the particular linear combinations 
of frequencies and momenta involved. Introducing a symbolic notation where 
uj denotes a positive linear combination of external frequencies and q a linear 
combination of momenta, the 2n-point functions have the scaling form 



v'p [ijj + ivpq)'^^^'^^ ' 

where an average over hot spots is understood, which leads to a real positive 
function of the frequencies and momenta represented by u and q. For the 
purpose of scaling, we have substituted J ~ ^/vp and 7 ~ g'^/vp. 

To estimate the relevance of the vertices in the scaling limit cu"^^^ ~ — >^ 0, 
where q dominates u in the denominator for z > 1, consider the 0^" term in 
the effective action [7]: 

where g2n is the coupling strength related to the vertex function b2n- Using the 
scaling dimension of the field [0^] = —{d + z + 2) (in frequency and momentum 
space) in two dimensions, 

[g2n] = -{2n ~l){2 + z)-[z- 2{n - 1)] - n{-4 - z) = {2 - z)n . 

The scaling dimension of all 2n-point functions is zero for z = 2, i.e., all 
bosonic vertices are marginal in the scaling limit, and it is not clear how to 
perform calculations with such an action. Our aim is to see if and how the 
disorder present in real systems changes the scaling dimension of the fermionic 
loops. 



2 Disorder corrections to a single fermion loop 

We consider static impurities modeled by a random local potential with mean 
squared amplitude [lU] ■ As long as no spin- vertex insertions appear between 
impurity scatterings, the disorder corrections in the Born approximation have 
the standard form. In the fermionic propagator G{ie,p) = {H — ^p)~^ the 
Matsubara frequency e is cut off as e = e + sgn(e)/(2r) at the scale of the 
impurity scattering rate 1/r = 2t:pqv?. Although the particle-hole bubble 
B{ie + iuj,ie,q) with small momentum transfer q is cut off by disorder, the 
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direct ladder resummation has the diffusive form L{ie + iuj,ie,q) ^ v? / {\uj\t + 
Dq^r) if both frequencies he on different sides of the branch cut on the real 
line. The diffusion constant is D = VpT/2 = VFi/2 in d = 2. Throughout this 
work we assume that impurity scattering is weak, l/ikpd^) = l/(ei?r) -C 1. 



2. 1 2-point function 



Due to the linearized dispersion around the hot spots, the bosonic self-energy 
is unchanged by disorder corrections to the fermionic propagators: 




spin - 27r (27r)2 ti - i(e + Uj) - ^p+g.+q 

= —7 \uj\ . 

As the direct impurity ladder with large momentum transfer is not singular, 
the leading vertex correction is given by a single impurity line across the 
bubble: 



X u 



1 

-7- , 

r 



27r (27r)2 ie - ie - ^p+q^ 
d^p' 1 1 



(27r)2 ie -ip ie- ^p/+q^ 



where the cutoff scale ei? is used for the divergent frequency integral. The mass 
correction shifts the position of the critical point as a function of the control 
parameter by a finite amount. We assume that the system can be fine-tuned 
again to the critical point by adjusting the control parameter. 

The correction to the bosonic self-energy from the maximally crossed ladder 
(cf. equation ( IB. II) ) is 

Sdirty V^^.Qc + Q)^1M • 

The logarithm does not change the bare scaling dimension, however the ab- 
sence of a corresponding term in q could tend to increase z and possibly make 
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the clean vertices irrelevant starting from z = 2. We will not further discuss 
this possibility and eliminate the crossed ladders by a small magnetic field 
because, as we shall see, the vertices in the presence of impurities will acquire 
new singularities due to diffusive ladders which will not be affected by the 
value of z. 



2.2 ^-point function 



Including only self-energy disorder corrections cuts off the frequencies in the 
denominator by cD = a; + sgn(ct;) / r and yields the scaling form (where u and q 
represent the same linear combinations of external frequencies and momenta 
as in equation ([T])) 



JO) 

^^4, dirty 



-nJg 



a 



UJ 



lUJ 



9_ 



l^tu + 1/r + ivFq)^ 



For vpq 1/t., this contribution vanishes linearly in uj and is, therefore, 
irrelevant in the scaling limit. Instead, a singular contribution is obtained by 
including diffusive ladders into the loop. Since direct ladders between two 
propagators with momenta separated by are not diffusive, we insert direct 
ladders only between propagators with almost equal momenta — i.e., near the 
same hot spot (see Fig. [3]). 




Fig. 3. The four-point vertex with one direct ladder insertion factorizes into pairs 
of spin vertices with three propagators between them [R vertices) and the direct 
ladder (L). 

The Raaa subdiagram is constructed from two propagators near one hot spot 
a and one propagator near the associated hot spot a, with two spin vertices 
in between: 



1 



1 



^ 9 0(— e[e + uj\) , ,, 



*e - G 



(3) 



which saturates to a constant for small tu, q. Combining the parts, we obtain 
for the 4-point vertex with one direct ladder (by the superscript we denote 
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the number of ladders) 



^Sirty = -J2 J ^ ^(^1 + 92, g + qi) L{qi + ga) R{qi + q2,q + 92) 

spin 

^^-^ (4) 

as the scahng form for typical external uj, q. We have checked explicitly that 
symmetrization of the spin insertions does not change the leading singularity. 
Note that the last factor in the ladder contribution (jl]) becomes larger than 
unity only for ur^Dq^r < l/lkpi), while ^i^dirty is cut off for ujT^Dq^r < 1. 
This appearance of two crossover scales is a central observation of our work. 
Adding the contributions with zero and one ladder, 

, 9\\iZH^ (^r and/or g£ » 1) 



^4,dirty ~ ^ ^ S 9 



(uT and o£ <C 11 



There are, of course, many more ways to insert ladders into the 4-point vertex. 
Since we are interested in the most singular contribution that determines the 
scaling, we have already excluded direct ladders between different hot spots 
on the basis that they are not diffusive. However, one could also add a second 
ladder in Fig. [3] between the a and f3 lines (with a = (3). This does not change 
the singularity, but gives a much smaller pref actor 1 / {kpiy. In general, all such 
crossings of direct ladders yield internal integrations over the ladder momenta 
(as in the case of the cooperon ladder) and give at most logarithmic corrections 
but not a higher singularity. Logarithmic corrections do not change the scaling 
dimension, and because of the higher order in l/{kpi), all diagrams with 
crossings of direct ladders will be neglected. Note that the numerical values 
of the coefficients of course depend on these diagrams. We argue further that 
the insertion of crossed ladders leads at most to the same singularity as those 
with direct ladders, up to logarithmic termsOl 

In conclusion, this leaves us with a much smaller set of diagrams displaying the 
leading singularity: all possible insertions of direct ladders between propagators 
of the same spin (at the same hot spot), which do not cross each other. 

For the 4-point vertex, the ladder diagram in Fig. [3] (with summation over hot 
spots and symmetrization of external lines understood) is the only one meeting 
these conditions and is therefore sufficient to obtain the leading singularity. In 
the scaling limit ~ g^, the ladder correction has the same scaling behavior 
(constant) as the clean vertex g^, hence the 4-point coupling remains asymp- 



^ Alternatively, as discussed above in connection with S^dirty^*^ '-'^^ could apply a 
small magnetic field, which does not cut off the singularity in the clean case and 
with direct ladders, but suppresses the contribution from crossed ladders. 



9 



totically marginal even when disorder is included. This leads to the question 
how the higher 2n-point vertices behave. 



2.3 Higher bosonic vertices 




Fig. 4. The disordered 2n loop (shown here for n = 3) with n — 1 ladder insertions 
arranged in a chain. This can be extended for larger n by repeating the (XL) block. 



Following the above discussion, we consider a particular insertion of ladders 
into the bare 2n-point vertex which meets the above condition for a maximally 
singular contribution. As shown in Fig. HI we group two adjacent spin vertices 
together (the R part as in Fig. [3]), 



R ^ t T 



RL ~ « 



This is connected via a direct ladder L to an X vertex made from four fermion 
propagators with two spin and two ladder insertions, 



X ~ — (for small external cu, q) 

Vp 



There is another contribution to X not depicted in Fig. HI with two spin 
insertions on the same fermion line: this term is of the same order of magnitude 
but generically does not cancel the one shown. 

Repeating the {XL) part n — 2 times and finishing with another R vertex, we 
obtain a chain-like diagram with n — 1 ladders which, for small u and q, can 
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be estimated by the scaling form 



■^2n, dirty 



spin 



9' 



n-2 o 

I 9 



9_ 



u + Dg2) (e^) + Dq^) (ep) J 

2n 



UJ 



Diagrams with fewer than n — 1 ladders have a weaker singularity and give an 
irrelevant contribution to the coupling in the scaling limit. 




L 



H„ 



V y 




Fig. 5. The disordered 2n loop for even n (shown here for n = 4) with n ladder 
insertions and a Hikami vertex Hn in the middle. 

There are other diagrams with n ladders, which at first appear to be even more 
divergent but upon closer inspection turn out to have the same singularity. 
As shown in Fig. [5|, for even n one can connect n R parts via n ladders to an 
ra-point Hikami vertex (TT| . 



^ r" (c. + Dq') 



The additional factor of an inverse diffusion propagator in Hn is due to the 
insertion of further single impurity lines which cancel the constant term and 
leave only terms linear in u and q^ for small u and q; this effectively cancels 
one of the n diffusive ladders. Connecting the RL parts to the Hikami vertex 
and adding the frequency integration along this one large fermion loop, 



'^2n,dirty 



X: / deHniRL)" 

spin ' 
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2n 



UJ 
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has the same singularity as the chain-type diagram. 

In appendix [Xj we propose a non-hnear a model for the spin modes coupled 
to low-energy diffusion modes (instead of the original electrons) with only one 
local (constant) coupling. This allows us to control all single-loop diagrams 
exhibiting the leading singularity, thereby supporting the perturbative calcu- 
lations in this work. Elimination of the diffusive modes would lead again to an 
effective action for the spin modes with infinitely many marginal couplings. 



3 Scaling of the disordered single loops 

3.1 Smallness of ladder corrections and second crossover scale 

As for 64,dirty5 the disordered loops (beyond the 2-point function) feature 
two crossover scales, one where disorder corrections in the self-energy of the 
fermion lines cut off the vertices, and another where ladder corrections lead 
again to marginal scaling. The existence of these two scales can be traced back 
to the presence of hot spots. 

For comparison, consider the case of forward-scattering bosonic vertices. Self- 
energy disorder corrections become important and cut off the fermionic prop- 
agators at ujT 1 and g£ ~ 1. Adding one disorder ladder implies adding also 
two fermionic propagators with nearby momenta and performing one momen- 
tum integration. This additional contribution can be estimated as 

u^^{l + u:Ty + {qlY r d'p ^{l + cory + iqiy 

^/(l + ury + {qiy - 1 J (27r)^ y^(l + cur)2 + (g£)2 - 1 ' 

=27rpoT 

The ladder correction becomes dominant exactly at the same scale ut ^ 
qi ^ 1 where the self-energy corrections appear. Hence, there is only a single 
crossover scale between clean and dirty behavior. 

Also in the case of backscattering bosonic vertices, the self-energy corrections 
set in at ut ^ qi ^ 1. However, the ladder insertions are modified due to the 
presence of hot spots for the two additional fermionic propagators: 

. J y^(l + cur)2 + (g£)2 - 1 

In comparison with the forward-scattering case above, there is an additional 
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^ integration which effectively replaces one factor r by l/ep-, such that the 
diffusive ladders become dominant only at a second, lower crossover scale 
ujT, {qiy ~ l/{kpi). This means that the impurity ladder scattering of electron- 
hole pairs with nearby momenta is less effective by a factor of l/{kpi) if the 
particles are forced by spin-vertex insertions to be near hot spots between 
impurity ladders. 



3.1.1 Scaling regimes 

According to the above discussion, we can identify three different scaling 
regimes for the disordered 2?T,-loops: 

{UJT, (qif » 1) 
(1 » CUT, {qif » ^) 
» ..r, {qtf) 

Note that each pair of fermion-like propagators in the clean case is replaced 
by one diffusive propagator and an additional factor of e^^ in the disordered 
case, leaving marginal in the scaling limit. The contributions to the vertices 
which become dominant in the different scaling regimes are visualized in Fig. [S] 
below, and in this figure it is made explicit how the constraints on ujt and q^ 
are to be understood. 



02n,clean ~ (^+jt,^<j)2{n-i) 
-'2ra,dirty ~ \ "2n,dirty (i/^2)n-l 



(n-l) 

2n,dirty (w+Dg2)n-i(e^)i 



There are two different diffusive regimes in the model: a fast one for charge 
modes with diffusion constant D, and a much slower one for spin modes with 
Z//7 ~ D/{kpt). We can, therefore, look at two variants of the scaling limit 
oj ~ q^. For charge diffusion we have uj ~ Dq^ (see the red/dashed scaling line 
in Fig. [8]), and 

'^2n,clean ~ (^/^n-i (^T > 1) 
^2n,dirty ~ \ &2n,cutofr ~ (;^/^2')n-i (1 ^ CJT ^ ^) (5) 
&2n,ladder ~ (^,^'')n-i > ^r) 

shows a non-monotonous behavior. On the other hand, if the typical values of 
OJ and q are given by the spin propagators connected to the external legs of 
the fermion loops, we put 70; ~ z/g^ (see the blue/solid scaling line in Fig. [8]), 
and 

(^^ » fe^) 

^2n,dirty ~ ^ "2n,cutofr ~ (i/^2)„-l (^^ > U;r > J^^) (6) 

This scaling analysis suggests that even though the ladder corrections scale 
marginally for asymptotically low frequencies, there is a range of frequencies 
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and momenta where the singular clean vertices are already cut off and the 
ladder corrections are still small, such that the Hertz-Millis theory might apply 
in this zone. As we shall see in the following section, such a regime may be 
hidden by further contributions from multiple loops. 



4 Disorder corrections to multiple fermion loops 

The multi-loop disorder corrections arise from impurity lines connecting dif- 
ferent fermionic loops. In the simplest case, one takes n static particle-hole 
bubbles with two spin insertions each (mass terms) and connects them with 
single impurity lines. As the impurity lines do not carry frequency, there are 
only n independent frequencies in this 2n-point spin vertex. The correspond- 
ing coupling in the action has therefore a different scaling dimension than the 
single-loop contribution (with 2n — 1 independent frequencies) and is generally 
more relevant. In fact, the missing frequency integrations lead to a scaling as 
in the classical field theory and the Harris criterion ^12j applied to the bare 
model implies that such contributions are relevant in (i < 4. 

We first define the n-loop vertices in the disordered spin-fermion model 
and then discuss the energy scales where these additional vertices become 
quantitatively more important than the single- loop contributions. Let A[l^] 
denote the particle-hole bubble at arbitrary momentum transfer q (not just 
near and zero external frequency = in the presence of a particular 
configuration of the disorder potential V , 



where ^ is the hopping matrix, the spin-fermion vertex F = exp(— igx) is a 
diagonal matrix in real space (with spin indices suppressed), and the trace 
runs over spatial indices. We diagonalize ^ + V = ^k\^)^k{k\ and perform the 
e integration. 



We now specialize to the case = q^ and define the static connected 2n-point 



A[V] 
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vertex as 



^2n 




' disorder average, connected part 



We assume that for a generic large momentum transfer near and generic 
band dispersion, all have a finite and nonzero limit as a;,, |gj — q^\ — >• 0. 
The corresponding terms in the action 



have only n independent frequency integrations but 2n — 1 momentum inte- 
grals, stands for the running coupling with bare value ^2n- The constraint 
on the frequency integration leads to a scaling dimension 

[^2n] = —nz — {2n — l)d — n{—d — z — 2) 

= d-n{d-2) {z = 2) 

= 2 {d = 2) . 

Thus, infinitely many couplings 62n are all equally relevant. The singularity 
found in two dimensions is so strong that one expects that even in three 
dimensions A4 remains a relevant perturbation and destabilizes a Gaussian 
fixed point in the disordered model, which may be related to the difficulty 
encountered in reconciling experiments with the Hertz-Millis theory for d = 3, 

z = 2 [nm]. 






Fig. 6. Chain-type contribution to the vertex (shown here for n = 3) with zero 
momentum transfer on the impurity hues. 

In order to estimate the quantitative importance of the vertices with 
respect to the 62n vertices, one needs to find the contributions to A2n at 
lowest order in l/lkpi). For n = 1, the leading term is the mass correction in 
equation ([2]). For n > 2, leading contributions to A2n are given by chains of n 
A[V^]'s connected by n — 1 single impurity lines with zero momentum transfer 
(see Fig. [6]). The terminal bubbles of the chain (with one impurity line) are 
roughly g'^/vp, while the intermediate bubbles with two impurity lines are 
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approximately g'^/lvpep)- The complete chain can therefore be estimated as 



A 



Jin 



Vpep 



u 



,n-2 



in > 2) 



In the scaling limit dominated by spin diffusion, we can replace {q/kp)'^ by 
oj/ep, and the single-loop vertex functions scale as given in equation ([6]). 
Below a certain frequency scale cu, the relevant couplings will necessarily 
become larger than the marginal couplings g2n'-, in order to find this scale 
one has to compare 62n,ciean with A2n/to'"~^, to account for the missing n — 1 
frequency integrations in the couplings with respect to g2n'- 
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Hence, in the spin scaling limit the higher vertices A2n start to dominate 
the clean single-loop vertices 62n,ciean at successively lower frequency scales. 
Likewise, the cutoff vertex 62n,cutoff is of the same magnitude as A2„/co'"'~^ for 
cur ^ l/(kp£)^'^"'~^^^^. Extending these arguments beyond the scaling limit 
to the uj-q'^ plane, we obtain the crossover lines between single- and multi- 
loop contributions indicated in Fig. [HI In summary, A4 dominates the clean 
vertex 64 below u < vpq/ ^fkpl for vpq >1/t, and the cutoff vertex 64 below 
^ ^ ^/{.T^kpt) for Vpq < l/r. The higher vertices A2n>4 become dominant 
for Vpq < 1/r only below oj ^ l/{Tkpt). 





(a) 

Fig. 7. (a) Ring-type contribution to the A2n vertex (shown here for n = 4) which 
is peaked when the impurity lines carry a momentum near 2kp. (b) Points on 
Fermi surface connected in the singular q^-2kp bubbles, with a commensurate 
and incommensurate 2kp. 
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In addition to the chain-type diagrams above, there are diagrams with n bub- 
bles arranged in a ring and connected by single impurity lines. In this case, one 
has to integrate over the momentum carried by the impurity lines, involving 
also momenta near 2k p where the bubble with two static spin insertions 
and two static 2k p charge insertions becomes singular in the clean case (see 
Fig- El)- This singularity is related to the well-known 2k p singularity of the 
particle- hole bubble [TS]; while in the clean case only a small region of the 
Fermi surface around the hot spots is visited, impurity scattering visits the 
whole Fermi surface, including the parts separated by 2kp. The disorder cor- 
rection to the fermionic self-energy provides a cutoff for this singularity, thus 
changing the estimated power of l/{kpi) in the expression for A2n- We esti- 
mate that all such diagrams are at least of order l/{kpiy, which implies that 
they become dominant at the same scale as Ag, i.e., only below u ^ l/{Tkpi). 

Therefore, in the spin scaling limit we still obtain clean anomalous behavior 
above ~ 1/r, which gets modified by a single relevant vertex A4 for 1/r > 
uj > l/i^rkpi). For u < l/^rkpi), ever higher multi-loop vertices A2„, add to 
the single-loop vertices. 
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Fig. 8. (Color online) Regions in the co-q^ plane where different contributions to the 
spin vertices become dominant. The lines separating the different regions are meant 
as a guide only, as they are determined only up to prefactors of order 0(1). 

The problem of a disordered AFM was considered in Ref. [16], where the au- 
thors analyzed an AFM 0^ model with a local interaction and a random 
mass term within an e = 4 — d expansion. Averaging over random mass 
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configurations yields a A(0^)tj(0^)a;' term from tlie static mass-mass corre- 
lation {mm), which depends on three independent momenta but only two 
frequencies. Therefore, A is a relevant coupling in c? < 4 and the clean model 
is unstable against arbitrarily small disorder. Recently, this model has also 
been studied using a strong-disorder RG [H]. 

However, according to the above discussion the clean AFM 0^ model is not 
applicable in d = 2 because u is nonlocal and there are infinitely many ad- 
ditional marginal couplings which even change the direction of the RG flow 
for M f7]. Furthermore, we have pointed out that for a generic fermionic band 
dispersion the average over a random fermionic potential implies a fluctuat- 
ing mass term which has non- vanishing higher cumulants (m"). While these 
higher cumulants are irrelevant in c? = 4, in two dimensions all couplings A2„ 
are equally relevant and may possibly change the direction of the RG flow also 
for A = A4. Thus, it is not obvious which of the results of [16] hold m d = 2. 



5 Discussion and conclusions 

We have shown that the fermionic loops with large momentum transfer, which 
are relevant to describe an AFM transition, exhibit two crossover energy scales 
if disorder corrections are added. In the scaling limit set by charge diffusion 
u ~ Dq^i the singularities of the clean (marginal) vertices g2n are cut off at 
scale a; ~ 1/r and the 4-point vertex and beyond become irrelevant. Below 
uj ^ l/(r/cp'£), however, diffusive ladders lead again to marginal scaling, albeit 
with a diffusive functional form of the vertex functions different from the clean 
case. On the other hand, in the scaling limit set by the slower spin diffusion 
oj ~ ^q^ ~ Dg'^/^kpi), these crossovers occur at the same momentum scales 
but at frequency scales smaller by a factor of l/{kpi). A summary of crossover 
scales and boundaries is reported in Fig. [HI 

In addition to the single fermion loops, there are diagrams made of multiple 
fermion loops connected only by impurity lines, which were not present in 
the clean case. In accordance with the Harris criterion applied to the bare 
model, these are relevant perturbations in d = 2 which make the clean model 
unstable against disorder. Indeed, in contrast to d = 4, in two dimensions the 
bosonic model contains infinitely many equally relevant vertices A2n due to 
disorder corrections to multiple fermion loops, which we estimate to become 
important below a; ~ l/^ry/kpi) (charge diffusion) or a; ~ 1/r (spin diffusion), 
respectively. 

Combining these results, we can distinguish two cases: if scaling is dominated 
by charge diffusion, the anomalous clean behavior |6| above u ~ 1/r is cut 
off below to make place for an essentially non-interacting behavior (Hertz- 
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Millis theory with only irrelevant couplings) until uj l/(r-\/Fp7), where the 
relevant disorder vertex A4 starts to dominate. Below uj ^ l/lrkpi), infinitely 
many more disorder vertices dominate the respective single-loop vertices. 

In the case of scaling determined by spin diffusion, the clean anomalous behav- 
ior is modified below ^ 1/r by a single disorder vertex A4, while the higher 
single-loop vertices 5'2n>4 are not yet strongly modified by disorder. Only below 
UJ ^ l/{Tkpi) the single loops are cut off, but at the same time more disorder 
vertices A2n>4 appear. In consequence, there is a direct crossover from clean 
anomalous to strongly disordered behavior which implies that the theory of 
Hertz and Millis may not be applicable for the two-dimensional antiferromag- 
net. 

TE wishes to thank A. Chubukov, A. Rosch and M. Salmhofer for fruitful dis- 
cussions. We thank the Alexander von Humboldt foundation (TE and CDC), 
and the Italian Ministero dell'llniversita e della Ricerca (PRIN 2005, prot. 
2005022492) for financial support. 



A Non-linear sigma model 

In this appendix we outline the steps to formulate the problem of disordered 
single loops in terms of the non-linear a model for interacting disordered 
electrons [T7|18|19ll2U] in order to support the perturbative calculations in 
this work and control all contributions with the leading singularity. Assuming 
some familiarity with the non-linear a model itself, we introduce only the 
modifications necessary to accommodate the spin-vertex insertions. We start 
with noninteracting electrons in the presence of disorder, add the spin vertices 
as couplings to an external field, perform the disorder average using the replica 
method, integrate over the fermionic degrees of freedom and obtain an effective 
action for the Q matrices (in standard notation), 

S[Q, 0] ^ I {fa Tr Q2 _ 1 Tr ln(Go ^ + - <p)} 

where Q are matrices in frequency and replica space with a weak real-space 
dependence representing the electron-hole pairs while 0(r) = 0(r) exp(irq^) is 
the staggered external field, with the matrix slowly varying in spaceLH The 
Q matrices can be expressed as a rotation Q = T'^QspT of the saddle-point 
solution Qsp = sgn(e) of the classical action for = 0. In the vicinity of the 
saddle point, one obtains 

S[Q,(i)] ~ J dr[DTT{VQy - 4Tr(eQ) - (terms in and Q)} , 

^ For simplicity, we do not write explicitly the frequency dependence of 0. 
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where the first two terms are the standard non-hnear a model for disordered 
electrons and additional terms are obtained by expanding the logarithm in 0. 
One thus obtains vertices {QipY which couple the diffusons to the external 
spin field. 

We parametrize the Q matrices asQ = e^/^QspC"^''^, where the diffuson prop- 
agator {WW) is represented by the direct ladder L. There is a term QspCpWcj) 
in the action which corresponds to the R vertex in equation ([3]), and terms 
W(f)W(f) and W^cj)'^ (corresponding to the X vertex) with two diffusons con- 
nected to two spin- vertex insertions. Higher vertices {Q4>)'^ with more than 
two spin insertions generate sub-leading contributions by the scaling argu- 
ments presented below. 

For the leading singularity it suffices to consider all possible ways to connect X 
and R vertices via diffusons, with the possible inclusion of the Hikami vertices 
Hk [H], which represent the interaction of the diffusons in the absence of 0. 
Each vertex contributes a 6 function of all momenta, while each ladder implies 
an integration over its momentum. Hence, the power counting depends only 
on the number ns of additional 6 functions beyond the overall 6 function of 
external momenta, which is the number of vertices minus one, ns = ny — 1 = 
nR + nx — 1. The Hikami vertices do not contribute to ns since they scale as an 
inverse ladder. The number 2n of spin insertions determines n^+nx = n, such 
that [g2n] = z — 2ns = z — 2{n — l) in accordance with our previous calculation. 
This maximal scaling dimension is valid for a large class of diagrams, two 
representatives of which are the chain-type diagram in Fig. H] and the star- 
shaped diagram in Fig. O 

The class of diagrams with leading singularity contains, however, contribu- 
tions with a different relative importance measured in powers of l/lkpi). The 
dominant terms of correspond to a single fermion line and no 

crossings of direct ladders, otherwise an additional factor l/{kp€)^ is generated 
according to the formula 

fe>i ^ 

which is valid for any connected diagram with ri/^^fc ^ Hikami vertices 
and < riij < vertices R as well as nx = n — n^ vertices X. For instance. 
Fig. [3] corresponds to RLR and has I = (with uh^^ = and riR = 2), while 
an additional ladder crossing the first one corresponds to Tt{XLXL) with a 
single closed fermion loop which has / = 1 (with = 0). 

If we include the spin propagator in the action and integrate over the dif- 
fusons, the purely spin-wave action is recovered with its infinite number of 
marginal vertices, while integrating over will likely lead to an action with 
an infinite number of marginal diffuson vertices. It appears that the action 
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written in terms of both (j) and diffusons provides the simplest formulation 
of the interaction of the low-energy spin and diffusive modes, with only one 
(constant) coupling associated to the coupling term / dr Tr [(50(50]. 



B Crossed impurity diagrams 



In this appendix we consider the disorder corrections to the bosonic self-energy 
beyond the Born approximation due to maximally crossed ladders Lc- In con- 
trast to the direct ladders they can have a diffusive contribution also between 
two propagators separated by an incommensurate (e.g., in the bosonic self- 
energy), and as the ladder momentum is integrated over, this typically results 
in a logarithm: 

X L^iu, p + p') G7+^(e + u) G7(e) 
(a;r«i) ff' 1 / 4 -ln(|a;|T) \ 

vj^ kpi ' ' [Ntt 1 + Dq^r J ' ^ ' ' 

While the above expression oc |a;| ln(|a;| r) vanishes for |ti;| — 0, it is logarith- 
mically larger than the Landau damping term oc and has the same sign. 
Thus, Ecrossed cuhauces the frequency-dependent part of the bosonic propaga- 
tor while leaving the momentum-dependent part unchanged for small frequen- 
cies and momenta. Therefore, there may be a tendency to increase z beyond 
z = 2. On the other hand, for N = 8 the term in parentheses becomes of 0(1), 
i.e., comparable with the direct ladder contribution, only if |a;| < lO^^r. As 
previously mentioned, in this work we discarded these contributions assuming 
that a small magnetic field would cut off this logarithmic singularity. 
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